function x=findgammass(gam)
global lamda Delta1 Delta2 sigma c1ss betta rss taumax Fkkss

muss=lamda*((1+(Delta1+Delta2/lamda)*(1-sigma))*c1ss^(-sigma)-sigma*c1ss^(-sigma-1)*gam*(1-rss*(1-taumax)));

x=-muss + betta*muss.*rss - gam*c1ss^(-sigma)*betta*(1-taumax)*Fkkss;
 